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The formulation of the initial value problem for the Einstein equations is at the heart of obtaining 
interesting new solutions using numerical relativity and still very much under theoretical and applied 
scrutiny. We develop a specialised background geometry approach, for systems where there is non- 
trivial a priori knowledge about the spacetime under study. The background three-geometry and 
associated connection are used to express the ADM evolution equations in terms of physical non- 
linear deviations from that background. Expressing the equations in first order form leads naturally 
to a system closely linked to the Einstein-Christoffel system, introduced by Anderson and York, and 
sharing its hyperbolicity properties. We illustrate the drastic alteration of the source structure of 
the equations, and discuss why this is likely to be numerically advantageous. 



I. INTRODUCTION 



►~^ ! ^^ ongoing research effort, numerical relativity, attempts to solve the partial differential equations (PDEs) per- 

I I taining to the Einstein equations using numerical methods. A significant current focus is centered on the study of the 
ly-v relativistic two-body problem [|lj. The broad features of the dynamics are already approximately known, through well 

^^ , established semi-analytical approaches, including perturbation theory (see e.g., a review in [^) and post-Newtonian 
t^ ' approximations (see e.g., latest developments in Isl and references therein). It seems intuitively right that one should 
^D ', be able to leverage such knowledge and reduce the computational complexity of the problem. There have been recently 
a number of proposals on various non-standard physical approximations to the spacetimes involved in mergers [0-0 . 
A common theme among such proposals is the elimination of degrees of freedom which are thought of as not crucial 
O . for some phase of the computation. There is also presently activity on the mathematical notions underlying numerical 
^ '' relativity, with particular emphasis on the hyperbolicity properties of the equations (see e.g., |8|]) and the nature of the 
?H . propagation of the constraints. Even while the mathematical structure of the general problem is under scrutiny, there 
PP' are particular instances of initial value problems which provide significant insights. For example, characteristic evolu- 
tion of spacetimes that admit caustic-free foliations appear to be completely tractable p] . Si gnifi cant improvements 
have also been claimed using approaches which decompose the field variables in certain ways |jlO| , [ll| . 
. _ We would like to advance in this paper a different perspective for numerical relativity. We will avoid any physical 

i3 approximation to the equations, but we will try instead to pin down geometrically and analytically the system by 

" " introducing additional "guiding" structures in the spacetime manifold. The key concept will be a background geometry, 

constructed analytically, or semi-analytically, and encompassing the a priori knowledge about the system. The true 
dynamics will then unfold as a "correction" on the background. The proposal is promising on several grounds. It is 
easy to show that reformulating the evolution in terms of finite perturbations can lead to significant gains in accuracy 
for systems that never deviate significantly from the background. This is also borne out in preliminary computations. 
It is possible that the stability of computations will improve as well. This is much harder to ascertain concretely. 
We note that currently the main avenues for improving stability are sought in the direction of securing hyperbolicity 
properties and the convergence of the constraints, but this effort has not been conclusive. The relevance of source 
terms seems less well understood. 

The concept of a background geometry has already been used in analytic studies of initial value problems (see e.g., 
in |1^ ) , and it has also been claimed that the use of a background spacetime clarifies the task of defining energy and 
momentum quantities for the gravitational field p3| . The treatment presented in [ p2| is particularly enlightening for 
the motivation of this paper. There, two four-dimensional metrics are introduced on the spacetime manifold M . One 
of them, say g^^, is promoted to physical status by virtue of satisfying the Einstein equations. The other, g^^, is 
relegated to a background role. Its associated Christoffel symbols, T^^p, are used for the construction of a covariant 
derivative V^. It is this derivative with respect to which one then expresses all relevant differential equations, e.g., 
gauge conditions and the Einstein equations. The fundamental evolution variables are the inverse metric differences 



X'"^ — dig^"^ ~ 9^^) I where g is the determinant of the g^" metric. With the imposition of the harmonic gauge, 
V^X^'^ — , one obtains the reduced set of Einstein's equations, 

g'""^p^aX^'' + L.O.T.'''' = -2R'"' + g^'^R , 

where L.O.T.^'' denotes terms of order lower than two and barred quantities in the right-hand side are the Ricci 
curvatures and scalars associated with the four dimensional background metric. Two remarks are worthwhile here. 
First, if g^i, is any exact solution of R^'^ — then the RHS of this equation is zero, which reflects the importance of 
the background for the structure of the source terms. Second, the hyperbolic characteristic structure of the equation 
is transparent, given its resemblance to a system of wave equations. It is noteworthy that whereas the source structure 
is background dependent, the hyperbolic structure solely relies on the physical metric. 

Following this example, we will formulate here the Einstein evolution problem with a view towards numerical 
applications. To this end we will work within the framework of the 3+1 decomposition. There are several reasons 
for opting for this path in view of the above mentioned (covariant) examples. Maybe most importantly, we would 
like to have explicit control over the gauge functions, as we would like to construct them using a priori information. 
Additionally, the 3+1 decomposition is tightly linked to a Hamiltonian framework. We will not explicitly exploit this 
feature here, but this will ultimately facilitate a systematic use of background spacetimes derived from post-Newtonian 
approximate solutions. 

The main result of this paper is the introduction of a set of new 3+1 evolution systems for the Einstein equations 
which explicitly separate a background geometry. The sets of evolution variables denote the finite deviations of the 
physical geometry from the assumed background. We derive systems involving second or first order spatial derivatives, 
with a corresponding change in the number of variables. In both cases we write the equations in a form which can 
explicitly cancel stationary or quasi-stationary balance terms, provided that adequate a priori description of those 
terms exists. We found that implementing the background geometry concept in conjunction with a hyperbolic system 
requires a slightly generalised geometrical construction. We will hence talk about the strong and weak versions of 
the proposal, according to the strength of the assumptions on the specifications of the physical and background 
spacetimes. All systems are constructed as three-covariant with respect to the background metric. 

The structure of the paper is as follows. In section II we review some results in order to establish the notation. 
In Section III we will describe a framework based on the introduction of a background spatial metric and develop 
the geometric objects emerging in this picture. We then write evolution equations for the spatial tensors describing 
the deviation of the physical geometry from a background one. The two natural sets emerging are closely related to 
the ADM and Einstein-Christoffel systems. We argue that the metric compatibility condition in the latter is more 
properly seen as a geometric identity on the three dimensional Riemann tensor. In Section IV we will generalise the 
framework slightly, in order to develop a symmetric hyperbolic system for appropriate evolution fields. In Section V 
we make a number of remarks. We discuss implementation issues and a number of potential applications. 

II. PRELIMINARIES AND NOTATION 

We will adopt the signature conventions used in |I4|. Throughout the paper, Greek indices denote spacetime indices 
running from to 3, whereas Latin indices are spatial indices and run from 1 to 3. Parentheses and square brackets will 
denote symmetrization and antisymmetrization respectively. Reference to geometric objects like vectors and tensors 
will refer to spatial objects, unless otherwise stated. We will systematically use an overbar to distinguish geometrical 
objects referring to a background geometry. Raising and lowering of indices will always be using the physical metric. 
We will consider only the vacuum case, i.e., absence of matter sources [T^^, — 0). Studies of matter systems (in 
particular fluids) along similar lines are pursued separately p5| and a generalisation will be given elsewhere. 

The 3+1 decomposition (also known as the Arnowitt-Deser-Misner (ADM) formulation) of Einstein's equations |lj] 
has been discussed in detail in many works (see, e.g., [0|). It assumes that the spacetime has topology i? x E, and 
hence it is foliated by the level surfaces S(t) of a scalar function t{x^). The unit normal to the hypersurfaces S(t), 
n^ = —NW^t, satisfies gfi^n^n'^ + 1 = 0. It can be used to decompose the four-metric as 

The vector field f^ (time threading of the manifold, t^^V ^t — 1) giving rise to the relations N = —t^^n^ and /3'^ = 
^M _ tVj^m allows to write the metric element in the form 

ds^ = -N^di^ + hij{dx' + I3'dt){dx^ + (3^dt) . 

Here, N is the lapse scalar, which is assumed to be positive, /?* is the spatial shift vector, and /i.y is the spatial metric 
(or first fundamental form) induced on the hypersurfaces S(i). 



In terms of the unit normal, the second fundamental form or extrinsic curvature of the spatial hypersurfaces is 
given by 

Kij=-^£nh,j => dthij ^ -2NKij , (1) 

where the time operator is dt = dt — £p, and £^ denotes Lie differentiation with respect to /?'. In particular, 

where (3i = hij(3^ . The Einstein equations are decomposed in a set of evolution equations for the components of the 
extrinsic curvature 

dtK.j = -D.DjN + N{R,j + KK,j - 2K,'K,i) , (2) 

the momentum constraint 

DjK',-D^K = 0, (3) 

and the Hamiltonian constraint 

R + K^ - K^^K,j =0. (4) 

In these equations K — h''^Kij, Di, Rij, and R are the canonical covariant derivative, Ricci tensor, and scalar 
curvature associated with the spatial metric hij respectively. The Einstein field equations are equivalent to the 
evolution equations (0) and (0) together with the constraints (0) and (0) . Given initial data on an initial hypersurface 
S(to), i-e. (/i° , K°A satisfying the constraints, and freely specifiable quantities N and /3* defined for all S(i), one can 
use the evolution equations to find the future development of such data. The framework is explicitly threc-covariant. 

III. 3+1 DECOMPOSITION WITH BACKGROUND METRIC: STRONG VERSION 

We assume a second spatial metric field, denoted by hij, cohabiting with hij on each of the hypersurfaces S(i). 
The four-dimensional manifold M is then endowed with an additional four-metric defined by 

where 

h^'''n^ = and h'^ = {h-^f . (5) 

The normal n^ = g^'^rii, is now a unit normal vector with respect to both metrics. Then, there is a unique projection 
operator onto the hypersurfaces S{i), h^^u = h^u = ^^y + n'^n^- As is clear, this construction uses the same lapse 
N and shift vector /3* for both spacetimes and hence the name strong version. Taking into account the relations 
N = —t^rif^ and /3^ = t^ — Nn^^ we write the background metric in the form 

ds^ = -N^dt^ + hij{dx' + P'dt){dx^ + P^dt) . 

The extrinsic curvature associated with the background geometry can be computed from the expression 

K^i, = -hi/Kf^anp , (6) 

where V^ is the connection associated with the background spacetime metric g^^- This leads to 

Kij = —^l^hij => Othij — —zNKij . 

The two spatial metrics give rise to compatible covariant derivatives Di and Di satisfying Dihjk — and Dihjk = 0. 
If r*jfc denotes the Christoffel symbols associated with hij then 8^ jk — ^''jk — ^^jk, the difference of the Christoffel 
symbols, is also a tensor. Indeed, under a coordinate change x* ^ x* = /* (a;*), the Christoffel symbols transform as 



p2 p 



dx'^ „, dx^ dx^ dx^ d 



a^i 



dx^ dxi' dx^' dx^ dx^' dx^' 



whence the tensorial transformation law for £^ jk follows 



C, jik' 



dx^ '' dx^' dx^' 



In fact, the explicitly covariant expression for f %fe is 



^*jfe — 2^* i^jkl + ^kjl 



where 



Hjk) , 



(7) 



Aki-i — Dhh 



^kij 



k '^ij 



(8) 



We see that A^ij is the covariant derivative of the physical spatial metric with respect to the background spatial 
canonical connection. The two covariant derivatives acting on a general tensor field are related by 



D,A'- 



D,A^-k...+£hiA'-' 



S'^kA'- 



I--- 



(9) 



Pictorially our geometric construct is presented in Fig. ^. The unit normal vector field n^^ represents Eulerian 
observers momentarily at rest on the hypersurface S(i) and is by construction the same for both spacetimes. The 
acceleration 1-form of those observers, a^j = n'^Vi/n^, (with a'^''nf^ = 0), is the same as seen from both spacetimes. 
Some calculations show that 

flp = /i^"V,(ln7V) = h^^'V^ilnN) = a^ , 

where we have used the fact that the background and physical four-covariant derivative have the same action on the 
scalar InA^. 



^t+dt 




FIG. 1. 3+1 decomposition and a background geometry. Two dimensions are suppressed for practical reasons. Two successive 
level surfaces of the time scalar are shown, with corresponding unit normal vectors in bold arrows. The standard decomposition 
of the time flow vector fleld into the lapse and shift components are shown. All those elements are independent of the intrinsic 
3-metric. The two different intrinsic geometries are also shown, as two different copies of hypersurface E(i) (which should 
be thought as overlapping). The filled circles on the hypersurface can be thought of as markers of equal proper length, thus 
indicating the different intrinsic curvature. The effects on the extrinsic curvature are also indicated. The right-most normal 
vector, when parallel transported to the left using the physical metric produces the dashed vector. The difference of this vector 
from the local normal (bold vector) is the physical extrinsic curvature. The same original normal parallel transported under 
the background metric, gives rise to the dotted vector, implying a different (background) extrinsic curvature. 



A. Time evolution equations and constraints 



Wc turn now to the task of writing the 3+1 decomposed Einstein equations in terms of a background metric. One 



further identity we will need relates the Ricci tensors of the two metrics (see, e.g. 
explicitly symmetric form of the Ricci tensor, such an identity is: 



12[) in terms of f*^.. Using the 



^ij ^ij 



= Dkr 



D{^£^j)k 



£^ 



j£ ki ^ £ iii£ 



l{i^ 3)k ■ 



(10) 



We rewrite the basic evolution system in the foUowing form: 

dth,j = -2NK,, , (11) 

+N£\,£'ki - iV£^(,£',)fe + N{KK,j - 2KMji) + NR,, , (12) 

where we have used the relation (|9|) to convert covariant derivatives with respect to the connection Di into covariant 
derivatives with respect to Di. Moreover, the Lie derivatives in the time operator dt are expanded as 

£f3Kj == P'^Dkh^j + h.kDjf]'' + hjkD^f3'' , 

and 

In the above equations, £^|j, is shorthand for equations (0,H). Hence this is the classic ADM form of the evolution 
system, but expressed using a background metric and connection. 

It will be useful here to introduce an alternative expression for Eq. (n3) . This form of the equation (for Gaussian 
coordinates) is given in |lq| and to write it we introduce a densitized extrinsic curvature (of weight 1) 

rj = hh'''Kkj, (13) 

where h is the determinant of the three-metric hij . The evolution of this quantity can be found through the definition 
of the extrinsic curvature (which implies the evolution equations dth — ~2hNK and dth^^ — 2NK^^) 

h-^hudtK}j = -NKK.,j + 2NK/Kji + dtK,, . 
Combining this expression with the evolution of Kij [Eq. (p^)] we get 

4/C^j = h^h''' {~DkD,N + Di{N£\j) ~ NDf^k£'j)i + N£'kj£"\rn - N £'\^^k£"' 3)1 +NRkj} , (14) 



Equation (h4J) departs from the standard form of ADM evolution systems, which explicitly preserve the symmetry 
of the extrinsic curvature and hence only require six components. Instead here one obtains the evolution to nine 
components, which satisfy the three additional constraints 

hi[ilC j] =0. 

A direct symmetrization, /C^ = hu^K} j\, e.g., after each evolution step in a numerical algorithm, would be sufficient 
for enforcing those constraints. 

An evolution equation for £L can be computed in the following way: first, using the evolution equation for hij 
[Eq. (|l])] and the Ricci identities for the commutation of spatial derivatives 

{D,D, - D,Di)A''-i... = i?^™,^"-/... + ■ • • - R^u.A'^-^,... - • ■ • , (15) 

where R'' jm denotes the Riemann tensor of hij , one can find that 

dtr.k = D^^Dk)p' - R\jk)il3' + D\NK,k) - 2D^j{NKk)') . (16) 

On the other hand, the Lie derivative of the Christoffel symbols is a tensor that can be expressed as follows ES] 

XpT'jk = DjDk(3' - R\,i[3' . (17) 

Combining Eqs. (|6|) and ( |l7| ) we get 

dtV'jk = D\NK,k) - 2D(^j{NKk)') . (18) 

which re-expressed for the background geometry reads 

dt£'jk = h'' {Di{NKjk) - 2^(j(iVi^fe),) + 2N£^.jkK^i] + ^f^fc , (19) 



where dtT^jk has to be replaced by the corresponding background version of Eq. (p^. 

With the promotion of the components of the tensor S^jk to independent variables, we must re-examine the addi- 
tional constraints those variables must satisfy. The first one concerns the metric compatibility of the connection, i.e., 
Dihjk — which, expressed in terms of the background connection, is equivalent to the definition oi E"^ jk [Eq. ([7|)], or 
to the following relationship 

Dkhij = 2f(ij)fe , where f^fc = huE^ jk ■ (20) 

Another condition, which can be seen as a consequence of (po|), is the metric- induced property of the Ricmann tensor 
of hij 

Rijkl + Rjikl — . 

In our framework, this condition expresses the fact that the action of the commutator of background covariant 
derivatives, [5^, Dj], when applied to the physical metric, is given by the usual rule 

(DkDi — DiDk)hij = —2hrn{iR"'' j)ki ■ 

n D[,l 
tensors (10[) in an explicitly symmetric form. 



This in particular contains the relation DuS'^jik = which has been used to write the differences between Ricci 



B. Non-linear perturbations 

In the previous subsection we introduced a number of systems which made explicit use of the existence of a 
background spatial metric. The spatial connection tensor S'^ij arose naturally with the consideration of a background 
metric. In contrast, we will need to explicitly introduce "perturbative" spatial metric and extrinsic curvature tensors. 
To simplify the notation we will use the symbol 6 to denote "perturbative" quantities: SA^'j... = A^'j... — A^'j.... 
Then, the relation between the "perturbative" spatial metric and extrinsic curvature tensors 

is given by 

On the other hand, we can always write the following evolution equation for the background extrinsic curvature 

dtK,, = -D.DjN + N{R,j + KK,, - 2K,'Kji) - S,j , (21) 

which has to be considered as the definition of Sij. In the special case where (h,K,N, (3) is an exact solution of the 
Einstein equations, Eq. (|2l| ) reduces to Sij = and can be interpreted as the dynamical evolution of an alternative 
initial data set (h°j,K°j). 

Subtracting Eq. (^l]) from Eq. (|2|) eliminates terms involving second derivatives of the lapse and the background 
Ricci curvature and introduces additional extrinsic curvature terms, 

^t^P^J = Dk{N£\j) ~ ND(^,£''j)k + N£\,j£\i - iVf ^(,f ',)fc + N8{KK,j - 2K,'Kji) + S,j . 



The extrinsic curvature terms can be absorbed with the use of the densitized variable lO j [see Eq. (^^] for both the 
physical and background spacetimes. Then, we can write the following alternative equation for the evolution of the 
components of the extrinsic curvature 

dt5V, = /ii/i''= [Di{N£',k) - iV(^(,£:'fc)z - £'jfc£"/„. + £'^(j£"'k)i) + midtK} , + S^k] , (22) 

where 

r]ij = -S{h^ih,j) = (h^i - h^i)hij - h^i<j)ij . 



It is worth noting that in equation (02) background Ricci terms have been eliminated. The choice of an exact 
background solution would eliminate Sij as well. The term involving dtJC"^j will be different from zero in general, but 
it is homogeneous in the metric deviations from the background. 



The expression of the evolution equation for E^ jk [Eq. ([l9|)] in terms of the "perturbative" and background quantities 
is the foUowing 



dtS' 



]k 









where ^j* — 6Kj^ and /i*-' must be understood in terms of (pij and hij, i.e. /i'-' = {[(p + h]~^y^ . 

In summary, utihsing the background geometry and identifying the key geometric variables we have derived four 
evolution systems, for the sets of fields {hij,Kij), (hij,K.ij), {(j)ij,ijjij,£''jk) and {4>ij , 51Cij , 6"^ jk) respectively. All 
those equations are formulated in terms of the background metric and connection, but the physical metric also 
appears explicitly. In addition, background source terms, e.g., Rij, dtt"^ jk are present. 

IV. 3+1 DECOMPOSITION WITH BACKGROUND METRIC: WEAK VERSION AND 

HYPERBOLICITY 



In what follows we introduce a modified version of the framework described in Sec. III. As before, in addition to the 
physical metric hij we assume a second metric field, denoted by ft,y , cohabiting with hij on each of the hypersurfaces 
S(t), but instead of considering the same unit normal vector we consider a new one n^ which we assume to be parallel 



to the physical one n^, i.e. 
the following way 



An^ (A > 0). We can combine these objects to define a background four- metric in 



hf"" 



l^ = hf'"' - y 



where /i^'^ satisfies the same conditions as in the strong version [see Eq. (|^)]. The condition that n'^ is a unit normal 
vector with respect to g^j^ = ft.^i/ — A~^n^n^ (hfi^n'^ = 0) reflects the fact that the lapse scalar for the background is 
different from that of the physical spacetime. We will call it N and assume is a strictly positive function. Then, we 
can take A = N/N and therefore the relation between normals can be written as 



N 



and 



N 

n' 



With the assumptions made in this weak version of the 3+1 framework there is still a unique projection operator onto 
the hypersurfaces S(i), h^i, — ft,^^ — S^\y + n^n^ — d^,y + n^^n^. 

Using the relations TV = —t^^n^ and (3^ — t^ — Nh^ — t^ — Nn^ we write the background line-element as 



o^ and (3^ 



^NMr 



Mj{dx'+(3'dt){dx^ +f3^dt). 

Pictorially, the current geometric construction is presented in Fig. 0. The Eulerian observers associated with 
the background metric follow the same worldlines as the corresponding observers in the physical spacetime, but 
their acceleration is different: Ui = Di{\nN) as opposed to a^ = _Di(lniV). Therefore, some of the rigidity in the 
correspondence between physical and background spacetimes is lost, but one still has a unique time evolution operator, 
dt, which is a consequence of having the same shift spatial vector /3' in both constructions. 




^t+dt 



FIG. 2. 3+1 decomposition on a background geometry: weak version. Two successive level surfaces of the time scalar are 
shown, with two sets of corresponding unit normal vectors (in bold and dashed arrows). The decomposition of the time flow 
vector field into the lapse and shift components is shown for the two possible lapse alternatives. The essence of the construction 
is that in order to advance from the level hypersurface E(t) to E(t + dt) one uses the same time operator dt ~ dt — £i3- 



As in Sec ._pi| , the second fundamental form associated with the background geometry can be constructed from the 
expression (|^7which leads to the usual relation between spatial metric and extrinsic curvature 



A. First order systems and hyperbolicity 



The evolution equations for {hij, Kij,£^jk), (|l l\) , ( [12D and (119), as they stand do not satisfy any hyperbolicity 
condition, which implies that we cannot formulate a well-posed initial- value problem (see MEOl). There are many ways 
of reformulating Einstein's equations in order to get a hyperbolic system of partial differential equations (see |p|,p3 25 



and references therein). The structure of the equations for (hij,Kij,S^ jk) suggests that the hyperbolicity properties 
will be very closely linked to the first-order symmetric hyperbolic formulation introduced recently in [E2|. This has 
been named the "Einstein-Christoffel" (EC) system, which is based in considering combinations of the variables 
{hij, KijjV^jk). In pq |, implementing this formulation by means of a pseudospectral collocation scheme, a stable 
evolution of spherically symmetric black holes with horizon excision has been achieved, and in ||27| , the total evolution 
time of 3D black hole computations has been extended by using EC type systems. 

As in the EC formulation [g^, instead of considering the lapse A^ as a freely specifiable scalar function, we will 
arbitrarily prescribe the scalar density (of weight —1) a, which in the case of the physical metric is defined through 
the following expression 

N{a, h) = ah^ , where h = det{hij) . (23) 

This automatically changes the principal part of the evolution equation for the extrinsic curvature [Eq. (|2|)] , which will 
be crucial for the hyperbolic structure of the final evolution equations (see [ESi for details). This way of prescribing 
the gauge precludes the possibility of choosing the background lapse N to be the same as the one of the physical 
spacetime A^, which is the main motivation for the weak formulation of the 3-1-1 framework given above. Instead, one 
can introduce a background densitized lapse a as above 

N{a,h)^ahK (24) 

We start the construction of the hyperbolic formulation with the evolution equation for the spatial metric "pertur- 
bation" (/)ij, which reads as 

40JJ = ~2N^,j - 25N K,j , where 5N ^ N ~ N , (25) 



and where A^ and A^ must be understood as given by the expressions (|23| ) and (24) respectively. Notice that now we 
have a source term proportional to the difference of lapses which did not appear in the previous formulation. 

With regard to the evolution equation for the extrinsic curvature "perturbation" ^ij , we must first rewrite the term 
DiDjN in term of a. Taking into account the expression for the covariant derivative of a scalar density a of weight 
w, DiG — diU — wT^ ija, it becomes 

DiDjN = 5(,(A^£^)fc) + h-2 {D,D,a + £\(,5,)a} . 
As is clear, this modifies the principal part of the equation for Kij (|lj). The result can be written as follows 

dtK.j = Dk{N£\j) ~ 2D^,{N£''j-)k) + N{£\j£'u + 2f [\(,f '1,,, + KK,, 

-2K/Kji + R,j) - h^D.Dja . (26) 

Then, the equation for V'ij takes the form 

4V'«j = Dk{N£\.j) - 2D^,{N£'^j-,k) + N(£\j£\i + 2£^\f^,£'\^)i) 

+5 [N{KK,j ~ 2KMji)] + dNR,j + h^D^DjU - h^D^Dja . (27) 

Finally, the equation for £'^ jk is 

dt^jk = h'' [DiiNK^k) - 2Df^j{NKk)i) + 2NK,^r^jk} 

+¥' {2D^^{NKk)i) - Di{NK,k)} ■ (28) 



The system of partial differential equations we have got for {(j)ij,ipij,£'^jk) [Eqs. (|2^), (p7|), and (p8|)] is still not 
explicitly symmetric hyperbolic. We can get such form by introducing a new quantity, with the same number of 
independent components as S^jk (eighteen), which will be considered as a fundamental variable |2|] 

Ckij — 2^kij — h^ih ™(A|;to|j) — Aj-ji^) — E[ij)k + hi^(ih ^{£\lm\j) — ^j)lm) ■ (29) 

The inverse relations are 

Afcij = 2 {C,kij + 'ihk[ih ™[Cj)im - C|im|i)]} , (30) 

£ijk = 2C(j7c)j - Cijk + '^hjkh ™C[ji]m • (31) 

Then, let us look for a system of evolution equations for the unknowns {(jjij^tpij, Qjk)- As is clear, Eq. ( [25| ) remains 
the evolution equation for (j)ij. The evolution equation for ijjij is obtained by replacing S^jk by Qjk in Eq. (Um. The 
procedure to carry out this calculation is to express first £'jfe in terms of Aijk [Eq. (||)] and then to use the following 
important property of Aijk 

D[iAk]ij — ~Kn(^iR™ j)ik , (32) 

which follows from the Ricci identities (|l5|). The quantity analogous to Aijk in the EC formulation is Qkij ~ dkhij- 
It satisfies the analogous simpler relation d[kQi]ij = 0. These relations constitute one of the keys to find a explicitly 
symmetric hyperbolic system. Once the property ( p2|) has been used, we have to express Aijk m terms of Cijfc 
[Eq. (po|)]. After some calculations we obtain 

— 2h' /l""' [CikmC,jln + 2{(kl{i — C{i\kl\){Cj)mn — C\mn\j)) — 2C(i|fem| (Cj)in " C\ln\j)) 
+ 2Cfe»,(0mn - Cmnl)]} + S [N{KK,j - 2K/Kji)] + S N R,^ 

+h^{£^^jDka - 2£^k{rDj)a - D.DjO) + h^D^D^a, (33) 



where £^ jk niust be understood as given by expression (31). 

To find the evolution equation for C,kij we have to apply dt to the expression ( |29|) and use all the information we 
have. It turns out that in order to have the desired form for the principal of this equation one needs to make use of 
the momentum constraint (H) that the physical spacetime has to satisfy. After some manipulations the result is 



^t(:k^J = ~Dk{NK,j) + 4iVi^,(,/i'"[g)i™ - C|/™|,)] + 27V/i,(ax,)'/i'»" - ,5,)'X™"]£,™„ 
+4/i^/ife(,<5,)'"/i'"i^i[,„[5„]a + a5^,]J + hk(,Dj0K) + hi(^Dj){NKk') 
+D\k\{Nkj)') - D\Nkj)k)] - /ife(,/ij>/i""[2A(iVXm") - W{Nki^)] . (34) 



At this point we have a system of first-order evolution equations for the unknowns {(f>ij , ipij , Ckij) [Eqs. (25), (|33| ) and 
(p4)], where some quantities constructed from {hij,Kij) appear explicitly but they must be considered as given by 
their expressions in terms of {cpij , -ipij) and (hij,Kij). We finish this section by showing the differential structure of 
the system ( p5| , p3| , p4[ ) . This analysis is standard (see, e.g., M) and one only needs to consider the principal part of 
the equations, which for our system looks as follows 

dtcfiij « 0, 

4V'y ~ -Nh'^'DkCuj , 

dtCkij ~ -NDk-ipij , 

where « means equality up to terms that do not contain derivatives of the unknowns {4'ij , ipij , Ckij ) ■ From the structure 
of the principal part we deduce that the system is symmetric hyperbolic. Indeed, the equations that determine the 
kernel of the principal symbol associated with the system are 

(n • ^)</«y = , 

(n • i)yj,, + i''Ck^J = , 

{n ■ OCkij + (.ki^ij = , 



where ^ = (^t,^i) is an arbitrary covector, n ■ ^ = g'^'^n^^„ = n'^^^, and ^' = /i^^Cj- From the structure of the 
principal symbol it follows that the characteristic directions of propagation of the system are the normal to the 
spacelike hypersurfaces, n ■ ^ = 0, and the light cone determined by the physical metric (n • ^)^ = 1, and therefore, as 
in the EC formulation, the characteristic speeds are the speed of light and the zero speed. From the analysis of the 
kernel of the principal symbol we deduce that the fields propagating along the normal to the hypersurfaces E(t) are 
all the components of 4>ij and 12 components of Ckij, namely 

Cfcij = ^fcSiSj" + ^-Ok{iS,j) + C'kij : 

where Ai, Bij, and Ckij satisfy: C^t = C^ij = CC^kij = 0. The fields propagating along the physical null cone are 
all the components of ipij and the 6 components of (^kij given by 

Ck^J = e{DikS.;S.j + Cfe^u + 26C(.^j)) , 

where now e = n • ^ (e = ±1) and the quantities Ei and Fij are such that ^'Sj — C^ij = 0- 

Finally, as it is usual in the theory of symmetric hyperbolic PDEs, we can introduce an energy norm for our system 
[Eqs. (11), dH) and (||)] which is given by 



E(i) 



S(t) 



b'^^,,+^'^^,j+C'''\Mj)d'S, 



where dS is the volume element of T,{t) and indices are raised using the inverse spatial physical metric h^^ . 

V. SUMMARY AND DISCUSSION 

The main results of this paper are the introduction of a set of new evolution systems for the Einstein equations which 
explicitly separate a background geometry. We have studied two broad classes, which differ slightly in their geometrical 
construction. An important realisation has been that for constructing a hyperbolic set of evolution equations one 
must relax the assumption of a common lapse. The reason lies with the fact that the use of a densitized lapse is 
important for inducing the appropriate modification to the principal part of the equations. This introduces a gauge 
function that incorporates information about the three metric. On the other hand, having a background spacetime 
which shares the same lapse and shift as the physical one is crucial for a maximal simplification of the resulting 
equations. The choices are summarised in Fig. y. The introduction of a background geometry makes all the systems 
explicitly three-covariant. Besides a certain aesthetic appeal, this feature may be handy in applications involving 
multiple coordinate systems, e.g., spherical coordinate patches adapted to black holes and infinity respectively. 




FIG. 3. The logical relationship between different possible arrangements for the background geometry. In the most general 
case (enclosing set), the background geometry has both its gauge functions (lapse and shift) different from the physical geometry. 
The use of a common lapse or shift vector defines the two inner sets. For the overlap of those two sets, both lapse and shift are 
the same as for the physical spacetime. This arrangement corresponds to the framework in Section II (strong version). The 
developments in Section III (weak version) refer to the set on the right, characterised by {N,l3,h). 
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We note that in determining the Hght-cone structure, or in the case of the first order systems, the propagation of 
eigenfields, it is the physical metric that is always the governing field. This implies that there would be no unpleasant 
side-effects associated with the choice of a background spacetime that is causally incompatible with the physical one. 
We have not elaborated on the Hamiltonian and momentum constraints. Those are essential for the construction 
of consistent initial data. The conformal techniques and decompositions involved in reducing those equations to 
well behaved elliptic systems seem at first glance not to be particularly benefiting from the assumption of a generic 
background geometry. 

A key issue that must be considered is the way in which the use of a background geometry would modify the 
numerical properties of evolution equations for the Einstein system. Common to most approaches to the evolution 
problem is the solution of a quasi-linear system of explicitly or implicitly hyperbolic equations with sources. For 
spacetimes involving compact objects, large curvature gradients are present. Those gradients survive in the stationary 
(non-dynamical limit) and involve a balance between flux terms in the principal part of the equation and the source 
terms. The source terms are not homogeneous in the evolution fields which can also be expressed as the existence 
of longitudinal modes. For fluid spacetimes (e.g., neutron stars) this balance has physical underpinnings in the fluid 
pressure supporting the star. For black holes it is geometric in nature. 

The implementation of the systems developed here can proceed in a number of manners. After a suitable choice 
of background geometry, a direct discretization of the equations using standard discrete approximations is possible. 
A practise usually adopted is to ignore the constraint equations in the construction of the algorithm (free evolution) 
and to monitor them as time dependent error norms. The integration of a physical data set will then consist of 
a sequence of discrete timesteps into the future. In principle, for a stable and convergent algorithm, the solution 
can be found at any future finite time to any desired accuracy given enough resolution. With the use of prescribed 
background geometry one has the flexibility of a more interactive approach to the time evolution, namely through 
an iteration of the evolution process. This technique could be useful in circumstances in which the background 
geometry is constructed analytically on the basis of a few poorly known parameters, e.g., post-Newtonian locations 
and momenta of black holes. The iterations would then provide improved estimates for those parameters. Whereas 
any of the iterative steps can in principle converge to the true solution, it has been our argument that the residual 
errors will be dramatically reduced as one improves the guess on the background spacetime. 

The incorporation of matter dynamics has not been discussed here. There is no specific obstruction to the process. 
In contrary, an implementation of closely related non-linear perturbation concepts to the relativistic Euler equations 
is already underway with promising results [iq ] . An optimal application candidate would be the study of non- linear 
perturbations of a single black hole spacetime. In this case the stationary exact black hole solutions provide readily 
available first guesses for the background geometry. More speculative, but well motivated, is the application to multiple 
black hole systems. Here the appropriate construction of an approximate spacetime would make heavy use of post- 
Newtonian solutions. The framework provides a well defined path for systematically incorporating post-Newtonian 
results into the full PDE evolution, which only hinges on the suitable expansion of post-Newtonian results into global 
metric data. 
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